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ABSTRACT 

An investigation 1s made of the statistical properties of blocks 
at the 500-mb level along the latitude belt 50-55N covering 180 de- 
grees of longitude including the block area. From the results ob- 
tained a method of computing a 72-hour forecast of block formation 
is developed. 

The author is deeply indebted to Professor Frank L. Martin 
for his advice and guidance during the progress of this investigation. 
Grateful acknowledgement is also due Professor Richard C. Camp- 
bell for his invaluable aid in programming the statistical routine used 


in the CDC-1604 computer. 
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le J Introduction 

Extended forecasting by statistical means is a field of meteor- 
ology relatively untapped yet showing vast promise of success. The 
most fruitful attempts appear to be 1n the field of weather-typing as 
used by Elliott!_ a et al, and weather typing by orthogonal functions 
as used by Malone [4 , chapter 28 | etal The one major shortcom- 
ing of weather-typing, however, lies in the difficulty of predicting 
the type to follow. 

The most persistent weather type in any One area is the block. 
Definitions of blocks vary widely but all definitions require a rela- 
tively stationary flow pattern for periods ranging from three to more 
than ten days. Predictibility of the formation of a block in any one 
sector would thus aid weather-typing considerably. Itis the purpose 
of this investigation, therefore, to devise a method of forecasting 
block-formations at 500 mb using statistical parameters. 

Rex 15] thoroughly investigated blocking patterns a the lower 
and middle troposphere, and correlated anomalies of temperature 
and precipitation with marked success. In his definition of a block, 
the flow pattern was required to maintain its identity for a period of 
at least ten days. Selfridge, Stevenson, and Wood [ 6] ie temutes 
referred to as SSW, based their findings on the definition used by Ser- 


ebreny ee that a block must remain in the sector of origin fora 





period of at least three days. In this study, all the characteristics 

of block types described by SSW and Serebreny have been utilized. 
Blocking-type boundaries are as defined by SSW and Serebreny. This 
investigation covers the period from six days before the formation of 
a block, F-6 day , until F-DAY, which is the day the block is actually 
within the defined type boundaries. 

Because of the availability of ocean stations and land stations both 
upwind and downwind, the Sub-Icelandic block was chosen for use as de- 
velopment data in this investigation. The formation area for the Sub- 
Icelandic block is contained within the area 50-65N and 10-30W. 

2. Data Collection 

Eleven stations were selected within the boundaries 50-55N and 
10E to 177W in order adequately to cover the area in which it was be- 
lieved indications of block formation would appear. The stations se- 


lected are as indicated in Fig. 1 below. 





Hol © ce. Weather plotting map illustrating stations and station indicators 


used in development. 
a 





For all fifteen winter cases of Sub-Icelandic blocks listed in 
the calendar of typesin SSW, the actual 500-mb data were collected 
for the period ranging from F-6 through F-DAY from the Historical 
Weather Map Series ion. 

The mean monthly contours from ''Normal 500-mb Charts for 
the Northern Hemisphere"'{ 11] were used to obtain the 500-mb mean 
height for each station in the month of each block. In some cases, 
linear extrapolation was used for computing the means. 

From the observed heights and their means, height anomalies 
were computed. These data were transferred to IBM punch cards. 
Linear correlation computations between height anomalies were then 
produced for all eleven stations taken two at a time, and for each day 


prior to and including F-3. 
3. Method of Computation 


All computations were made using the CDC-1604 computer. 


Data were programmed using the formulae listed below: 


2,  @. =3 (1) 
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In the set of formulae above, n=15 is the sample size, z is the 
De 
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(7) 


height anomaly at station x for the specified day before block formation. 
Le is a correlation between A and Zy- In most cases correlations for 
the same day were preferred (Table 1); correlations involving different 
time lags were also performed(Table 2). In addition, 24-hour height- 
change values were correlated across the entire network of stations 
(Table 3). In general, these gave many excellent significant correla- 
tions but were not considered to be as useful predictors as the previ- 
ous two types of correlation coefficients. 

In all computations, heights as measured, and anomalies as com- 
puted are in tens of feet. 
4. Evaluation of Computations 

From Hoel fi] , assuming that the long-term theoretical cor- 
relation of height anomalies between stations is equal to zero, and that 
the height anomalies are normally distributed, an actual value for the 
correlation coefficient (based on 15 independent cases) whose magni- 


tude 1s equal to or greater than 0.473 is significant at the 95% level of 





belief. The criteria for the selection of prediction stations were that, 
in all cases, the correlation coefficients between stations must be 
greater than 0.473 in magnitude and, furthermore, there must be at 
least three stations which are mutually and significantly correlated. 
Trial predictor stations were selected from the list of correlations 
shown in Table 1. 

Since it had been previously determined that a 72-hour forecast 
would be developed, the correlations for F-2, F-1, and F-DAY were 
not used in the development. 

Graphs of height anomalies were next drawn for each pair of 
significantly correlated stations associated with block formation. These 
graphs are shown in Figs. 3 through 9. 

Using a table of random numbers: le _ , fifteen cases of non- 
blocks were randomly selected and height anomalies for these cases 
were plotted on the appropriate graphs for the corresponding day. The 
regression line for block cases only was plotted on each graph. 

A discriminant analysis of the type used by Miller! 3 | was next 
performed on each graph. It was discovered, as expected, that the 
correlation between a station the the block area and one immediately 
adjacent was uniformly positive, and in no way discriminated between 
blocks and non-blocks. Of 17 graphs tested, 10 were rejected as hav- 
ing no forecast value. 

On the graphs which were retained, subjective analysis of the type 


used by Thompson I 9 was conducted. This consisted of drawing rel- 
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ative confidence isolines numbered from 0 to 5 on each graph. The 

"5'' isoline encircles clusters of block cases near the regression line, 
and the ''0"' line is near the axis of the perpendicular to the regression 
line. Essentially the areas within the '0" and ''5' lines test the hypo- 
theses to reject or accept block formation respectively, while the re- 
maining areas delineate the intermediate critical region No block case 
was permitted to correspond toa relative confidence value of less than 
three. This in effect assigned spacial weights in the critical region. 

The critical region was further tested in the following manner: 
(1) it is assumed that the events leading up to block formation are se- 
quential in nature and therefore to be equally weighted; and (2) the 
relative confidence value from each graph is summed for each day, 
thus assigning equal weight for each graph The graphs which display 
these relative confidence i1solines are shown in Figs. 3 through 9. 

It was noted that some blocking cases which had high relative 
confidence values on F-6 and F-5 had low values on F-3. Since no 
combination of pairs of stations gave good correlations on F-4, it was 
necessary to use a discriminant analysis between F-5 and F-4. Ac- 
cordingly, the relative confidence values for the four graphs on F-6 
and F-5 were summed and termed the first predictor (E,) . Since each 
day was to be equally weighted, the second predictor (E>) was formed 
by summing twice the relative confidence value of the F-4 graph plus 


the sum of the F-3 graphs. These two predictors for non-block and 
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Table 1. Table of significant height-anomaly correlations with no 


lag. 
Anomalies 
Correlated F-6 F-5 F-4 F-3 F-2 F-] F-DAY 
WA 0.47 0.68 0.56- Vt 7Oees OPS. OSS 
AB O66" 20455 0.60 0255.) -70e52 
BC 0.68 0.47 0252 0.53 
CD 0.49 
mo 0.49 
oR 0.79 0.69 O52 0.79 0.74 0.75 
RS 0709 0585 Nos O86, “0277 G52 SORSs 
o,f O55 0.49 
ea O255 
Vv 
WB 
AC 
iB D ~0.5] 
CO ~0.47 -0.63* -0.49 
DER 0.66 
Qs 0.58 0-56 
RT 
sU -0.55 -0.66 
coy ~0.62 -0.58 = 0577 
WC 


*x - Denotes those correlations used in the forecast procedure 
a 





Table 1 (continued) 


Anomalies 
Correlated F-6 F-5 fe oa F-3 


AD 

Be =O 20h 

Ck =O 39% 

DS 

Orr 

ee) -0.48 

SV =O257 
wD ~0.52 

AQ 

BR ~0.62 
CS 

mT ~0.49 

QU -0.49 

me V 

W OQ 

AR -0.58 


BS ~0.65 


QvV 


F-2 F-] F-DAY 


-0.47 


-0.49 


=0.65 -0.48 


=O 57 -0.49 


-0.62 


-0.50 


* - Denotes those correlations used 1n the forecast procedure 





Table 1 (continued) 


Anomalies 


Correlated F-6 


WR 


AS 


Bet 


Cru 


Dy 


Ws 


ee 


WU 


AV 


W V 


-0.70% 


-0.48%* 


F-5 F-4 Roo. eo 
-0 59% 
-0.49% 
-0.50 

-0.67 


F- | 


F-DAY 


-0.53 


-0.48 


™ — Denotes those correlations used in the forecast procedure. 





Table 2. Table of significant height-anomaly correlations with various 


lags. 

Anomalies Correlation 
Correlated Coefficient 
WoT, -0.54 
A3R, -0.58 
A358, -0.51 
SoT, 0.65 

2 -0.49 
C1Y~% -0.63 
A3D, -0.47 
WeAs 0.69 
AsD> -0.55 
C5Q,4 01-05" 
C.V, 0255 

D,D, O50 

DQ, 0.55 

RS. 0.49 

S.V> -0.71 

oe 0252 


= - Denotes those correlations used in the forecast procedure 
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Table 3. Table of significant correlations of 24-hour height 
tendencies with various lags. 


Anomalies Correlation 
Correlated Coefficient 
Bio Q54 0.90 

Bio C35 -0.50 

Can R54 -0.66 

Ws4 D423 -0.60 

Wey S39 -0.60 

Bo; ©32 Oc 

oO 243 053 

D3 Qe 4 0.54 

Qe 4 532 0.63 

Tey U25 =0751 
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block cases were next plotted (Fig. 10). The centers of gravity for 
block cases and non-block cases were next determined and a line join- 
ing these centers of gravity was constructed. The standard deviation 
for non-blocks was next determined and a distance along the line con- 
necting the two centers of gravity of 1.64 standard deviations from 
the center of gravity of non-blocks was marked. A line perpendicular 
to the line joining the centers of gravity was drawn through this point. 
Assuming the predictor values are normally distributed, non-blocks 
lying on the block side of this line have less than a 10% chance of oc- 
currence. 
5. The Forecast Procedure 

As stated earlier, it was assumed that the mechanisms or the 
series of events correlated are sequential in the case of the Sub-Ice- 
landic block. For the forecast procedure, it was further assumed that 
this series of sequential correlations holds true for any block and that 
only a shift in base line is necessary for computation. Therefore, a 
base line as close as possible to the one used in the Sub-Icelandic 
block computations has been developed. For simplicity, the latitude 
line of 50N was selected as the base line. Station locations are as in- 


dicated in Fig. 2 below. 


2 





DEVELOPMENT BASE LINE 





55.N 
T2038 -v2867 = 72848 
i e S 7206 | 
\ ? a eae 03322 
| l ; 72836 1 J @ 10338 
70454 1 , > ' Qg ° 
Q | | | | 
l 5" j 1 ] I | | | \ } 
| ' 50N 
| © | l 
pl l 27k ie e 
1 1 i . 4 | \ ail | 
; 1 ! | | | \ | 1 1 
V U be S R Q D Cc B A W 
2 54° sie 2159 eg 219) 9S 45° Selb So yg 17° 
~ * 
FORECAST BASE LINE 
V if S R Q G B A 50N 


Fig. 2. Comparison of spacing of station indicators used in both devel- 
opment phase and forecast phase. 


ior any type of block, the longitude of the center of the block- 
tye Area as defined by Serebreny [ 5,6! and SSW is used as the long- 
itude of station B. The other stations are located east or west of this 
point as indicated by the ''forecast base line" of Fig. 2. 

Using the same station identifiers as in the development stage 
but oriented relative to the new base line, height anomalies are compu- 
ted for stations A, B, and V for F-6 and stations C, Q, and R for F-5. 
Relative confidence values are determined from the appropriate graphs 
and summed for the first predictor, 

E, s c(Ag V¢) + c(B; V¢) + c(C,Qz) + c(CeRz). - (8) 
Preieht anomalies are then computed for station Q on F-4, and for sta- 
tions A, S, and T on F-3. Relative confidence values are thenobtained 


from the appropriate graphs and summed according to the formula be- 
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low to obtain the second predictor 
B= 2c(C,Q,) + c(A3S3) + c(A,T,). (9) 
Test data were thus obtained for each of 15 winter cases of blocks 
and non-blocks in the following areas: 


1) three cases each of Bering Sea-Western Alaska blocks, Sub-Aleu- 
tian blocks, and Scandanavian-Baltic blocks; 


and 


2) two cases each of Sub-Icelandic blocks, England-North Sea blocks, 
and Western Canada blocks. 


The results of these computations are indicated in Fig. 10. The 
data for these computations are compiled in Appendices III and IV. 
6. Results and Conclusions 

From the space correlations in SSW, the expected frequency of 
block days in all four sectors which they investigated is 903 block days 
out of a possible 2168 days. Of these block days, 209 were initial for- 
mation days. The remaining block days were rejected from the total, 
since only initial block days and non-block days were used in testing. 
The climatological expectency of block formations thus is 0.142 per 
day. The climatological expectency of non-block days is then consid- 
ered to be 1-0.142—0.858. From the sample summarized in Table 4, 
a forecast accuracy of 87 percent and a skill score, based on the cli- 


matological expectencies, of 0.73 were obtained. 
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Table 4. Contingency table of forecast results. 


i. observed 
ree) hice Mon hiocko 


eH ab a a EY 








Blocks - 

Non-blo cksi 2 oe ee 
In conclusion, although few cases were used to derive this fore- 

cast precedure, the results indicate that sequential correlations do ex- 
ist, both upwind and downwind during block formation, that can be uti- 
lized to make a 72-hour forecast. Itis believed that 1) by increasing 
the base line to a full 360 degree arc around the latitude belt of 50N, 
2) by decreasing the station-to-station interval to a constant 15 de- 
grees of longitude, and 3) by studying many more cases, that sequent- 
ial correlations will be found to exist that will correctly forecast not 
only block formation but also block dissipation. It1s also within the 
realm of possibility that the forecast can be extended to a period of 
greater than 72 hours. Itis further believed that such forecasting 
methods as presented herein may be of great value in determining em- 


pirical relationships to be entered as corrections or modifications to 


the results of numerical prognostic maps as currently being issued. 
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APPENDIX | 
Height anomalies for 15 cases of 


Sub-Icelandic blocks 


(Values given in tens of feet) 
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Case# 


10 


1] 


4 


13 


14 


15 


AG 
-25 
- 36 
-47 
=31 
+81 
+35 
-20 


-28 


-53 
+32 
+40 
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ceo) 


-18 
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APPENDIX II 
Height anomalies for 15 cases of 


non-blocks in Sub-Icelandic area 
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APPENDIX Ii! 
Test data; height anomalies for 15 block cases 


miscellaneous areas 


30 








Case Block Type 


els a Bering Sea-Western Alaska 
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36 Sub-Aleutian | 
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15. | Sub-Aleutian 
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APPENDIX IV 
Test data; height anomalies for 15 cases of 


non-blocks, miscellaneous areas 
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